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ABSTRACT 

A magnetic flux tube may be considered both as a separate body and as a confined field. As a field, 
it is affected both by the cyclonic convection (a-effect) and differential rotation (f2-effect). As a body, 
the tube experiences not only a buoyant force, but also a dynamic pressure due to downflows above the 
tube. When these two dynamic effects are incorporated into the aVL dynamo equations, we obtain a 
dynamo operating in the convection zone. We analyze and solve the extended dynamo equations in the 
linear approximation by using observed solar internal rotation and assuming a downflow suggested by 
numerical simulations of the solar convection zone. The results produce: (i) the 22-year cycle period; 
(ii) the extended butterfly diagram; (iii) the confinement of strong activity to low heliographic latitudes 
|$| < 35°; (iv) at low latitudes the radial field is in an approximately it phase lag compared to the toroidal 
field at the same latitude; (v) the poleward branch is in a 7r/2 phase lag with respect to the equatorward 
branch; (vi) most of the magnetic flux is present in a strongly intermittent form, concentraed into strong 
flux tubes; (vii) the magnetic field peaks at a depth of r = 0.96 Rq; (viii) total solar irradiance varies in 
phase with the solar cycle activity, having an amplitude of 0.1%; (ix) solar effective temperature varies 
in phase with the solar cycle activity, having an amplitude of 1.5 °C; and (x) solar radius also varies in 
phase with the solar cycle activity, having an amplitude of 20 mas. All these results are in agreement 
with the corresponding observations. 

Subject headings: Sun: interior — Sun: magnetic fields 
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1. INTRODUCTION 

Since Parker set down the foundations for solar afl dy- 
namo theory in his classic 1955 paper, many solar cy- 
cle models have been proposed, for example, the con- 
vection zone (CZ) dynamo models (Parker 1955), the 
overshoot layer (OL) dynamo models (e.g., Riidiger 
& Brandenburg 1995), the interface (IF) dynamo 
models (Parker 1975; Charbonneau & MacGregor 1997; 
Markiel & Thomas 1999), and the Babcock-Leigton (BL) 
dynamo models (Choudhuri et al. 1995; Durney 1995; 
Durney 1996; Durney 1997; Nandy & Choudhuri 2001). 
In CZ dynamos, differential rotation (the fJ-effect) in the 
solar interior shears weak poloidal magnetic fields into in- 
tense toroidal fields. Meanwhile twisting motions associ- 
ated with buoyancy, turbulent convection, and the Coriolis 
force transform toroidal fields back into poloidal fields (the 
a-effect). Since most of solar magnetic flux is present in 
flux tubes and they are subject to magnetic buoyancy in- 
stabilities, if there is no force to balance the buoyancy, they 
cannot be stored in the convection zone for a timescale 
comparable to the cycle period (Parker 1975). 

In order to circumvent the difficulty of flux storage in the 
main bulk of the convection zone, OL dynamos were pro- 
posed. Observations (Fig. 1) in fact show that most of the 
shear (the f2-effect) is concentrated in a thin layer near the 
bottom of the convection zone, known as the tachocline. 
This was seen as evidence that the tachocline is of key 
importance for the working of the solar dynamo. Such 
dynamos, use the same concept as CZ dynamos. In such 
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dynamos, the butterfly diagram comes out right, but the 
radial field is found to be nearly in phase with the toroidal 
field, instead of being in antiphase as observed, and the 
cycle period tends to be too short for thin layer models. 

Unlike both CZ and OL dynamos, in which the a- and 
O-cffects are assumed to be cospatial, in IF dynamos the 
a-effect operates on the high-diffusivity side of an inter- 
face, while the il-cffcct is limited to the low-diffusivity side. 
Their attractive feature is that the toroidal field generated 
on the low diffusivity side can be made arbitrarily strong 
by reducing the value of magnetic diffusivity there. The 
base of the convection zone, where the convection zone 
and the radiative zone are separated naturally, provides a 
physical interface. 

Unlike CZ, OL, and IF dynamos, which interpret the 
butterfly diagram as a dynamo wave, BL dynamos essen- 
tially work like a conveyor belt: the poleward meridional 
circulation near the surface transports the poloidal fields 
towards the poles at high latitudes, giving rise to the pole- 
ward branch of the butterfly diagram. At the poles, the 
fields are advected down to the bottom of the convection 
zone where the shear converts them into toroidal fields 
that get amplified while advected towards the equator. 
Once these are strong enough, they are supposed to form 
buoyantly emerging loops at low latitudes that give rise 
to active regions. The Coriolis force lends an inclination 
to the loop planes, introducing an a-effect to regenerate 
the poloidal field. In such models, the low-latitude con- 
finement of strong activity comes out naturally since the 



2 



toroidal field is amplified by shear as it is advected equator- 
wards, and the extended butterfly diagram is reproduced 
approximately in the sense that the latitude where the two 
branches part tends to be too high. A serious problem is 
that the approach only works with an unrcalistically low 
value for the turbulent diffusivity in the convection zone in 
order to keep the two parts of the conveyor belt separated. 

A successful solar cycle model should interpret the fol- 
lowing observational facts: 

1. The 22-year cycle period; 

2. The characteristic migration pattern (extended 
butterfly diagram, e.g., Makarov & Sivaraman 
1989); 

3. The confinement of strong activity to low 
heliographic latitudes |$| < 35°; 

4. The phase dilemmas (Stix 1976): 

(a) at low latitudes the radial field is in an 
approximately 7r phase lag compared to the 
toroidal field at the same latitude; 

(b) the poleward branch is in a 7r/2 phase lag 
with respect to the equatorward branch; 

5. Most of the magnetic flux is present in a strongly 
intermittent form, concentraed into strong flux 
tubes (Weiss 1964; Sheeley 1966; Stenflo 1973; 
Howard & Stenflo 1972). 

Helioseismic observations also impose rigorous constraints 
on the model: 

6. Solar internal differential rotation rate from 
inversions of helioseismic observations (e.g., Antia 
et al. 1998) should be used in the model; 

7. Helioseismic inversions made by Antia, Chitre, and 
Thompson (2000) indicate that there is a magnetic 
field of approximately 20 kG strength peaked at a 
depth of r = 0.96 Rq in 1996, rather than at the 
base of the convection zone. 

Furthermore, the model should reproduce the observed 
cyclic variations of total solar irradiance, solar effective 
temperature, and solar radius: 

8. Total solar irradiance varies in phase with the 
solar cycle activity, having an amplitude of 0.1% 
(Frohlich & Lean 1998); 

9. Solar effective temperature varies in phase with 
the solar cycle activity, having an amplitude of 
1.5°C (Gray & Livingston 1997); 

10. Solar radius also varies in phase with the solar 
cycle activity, having an amplitude of 20 mas 
(Emilio et al 2000). 

Using these constraints, we can see that none of the above 
solar dynamos are satisfactory. In particular, criterion 7 
rules out the last three models since they predict that so- 
lar magnetic fields peak at the base of the convection zone, 
leaving CZ dynamos as the unique candidate. 



Recently, Li and Sofia (2001) have shown that a mag- 
netic field of 20 — 50 kG located at r = 0.96i? Q can re- 
produce the observed cyclic variations of total solar irra- 
diance, solar effective temperature, and solar radius. If 
we use the helioseimic internal differential rotation in the 
CZ dynamo, the CZ dynamo will automatically satisfy 
criteria 6 — 10. In order to meet critcron 5, the mag- 
netic fields should be considered in the form of strong flux 
tubes. As a body, a flux tube in the convection zone ex- 
periences not only a magnetic buoyancy since gas in it is 
less dense than the surroundings (Parker 1975), but also a 
dynamic pressure if flow is present. Stratification provides 
a preferred direction; turbulence in the convection zone 
is likely to be anisotropic. In fact, numerical simulations 
of the solar outer convection zone (Chan and Sofia 1989; 
Kim & Chan 1998; Nordlund 1999) indicate a major pres- 
ence of downward-moving plumes. This suggests that the 
dynamic pressure of the downward-moving flow may push 
down the magnetic flux tubes to prevent them from rising, 
until they reach certain depth in the convection zone. 

In this paper, we extend the classical CZ dynamo model 
by taking into account the magnetic buoyancy and down- 
ward flow effects. We show in §2 that these two effects 
can be incorporated into the mean field dynamo equations 
by modifying the turbulent diffusivity (3. In §3 we lin- 
early analyze the dynamo equations to obtain the observ- 
able quantities of the model. In §4 we use the theoretical 
findings to interpret the observational facts listed above. 
Finally, in §5 we summarize the results. 

2. BASIC EQUATIONS 

We consider a magnetic flux tube with radius a and 
length L. Its magnetic field B is assumed to be uniform. 
As usual, we decompose it into a toroidal and a poloidal 
component, 

B = Be + V x (Ae ), (1) 

where is the azimuthal unit vector, B is the toroidal 
field component, and A is the (toroidal) vector potential 
of the poloidal field. The local hydrostatic equilibrium of 
the flux tube in the solar interior requires 

P = Pi + B 2 /8ir, (2) 

where P and Pi are the total and internal gas pressure, 
respectively. Assuming that this equilibrium in the flux 
tube is reached solely by adjusting its density, then the 
internal density pi is related to the normal density p by 

Pl = pPi/P. (3) 

As a result, the density reduction of the magnetic flux tube 
is 

p - Pi = pB 2 /8wP. (4) 

The reduced density inside the flux tube produces a buoy- 
ancy, 

F = ira 2 Lg(p- pi)e r , (5) 

where g is the gravitational acceleration. We assume that 
a is smaller than the pressure scale height so that we can 
consider that the gas density of the flux tube is uniform. In 
this case, the total mass contained in the flux tube equals 

Mi = na 2 Lp(l - B 2 /8nP). (6) 
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(9) 



Therefore, the buoyant acceleration g b can be expressed 
by the gravitational acceleration g, total pressure P and 
the magnetic field strength B as follows: 

g b = gB 2 /(8TrP-B 2 ). (7) 

If there is no downward flow to balance this pressure, 
the only force that goes agaist the buoyant force is the 
aerodynamic drag, 

F d = -\C d pu 2 aLe r , (8) 

where u is the velocity of rise and Cd(<~ 1) is the drag 
coefficient. Taking into account both the buoyant force 
and the drag force, as done by Parker (1975), the terminal 
velocity of rise occurs for F = Fd , yielding the rate of rise 

/ 2irag p- pi\ I ira 

U= {——) = Va {c^H p 
where H p = P/pg is the pressure scale height, Va = 
B/(4ttp)^ 2 is the Alfven speed. The rise time to the sur- 
face is smaller than the required field amplification time 
of about 10 years by the classical dynamo theory, as esti- 
mated by Parker (1975). 

However, numerical simulations of the Sun's outer con- 
vection zone made, e.g. by Nordlund (1999), indicate a 
major presence of downward-moving plumes of high ve- 
locity. This downward flow has an inward dynamic pres- 
sure pv 2 . If this pressure can balance the buoyant force 
of the flux tube, then the flux tube can be stored in the 
convection zone so that it can be amplified by the afi dy- 
namo. Antia et al (2000) employ the observed splittings 
of solar oscillation frequencies to separate the effects of 
interior solar rotation, and to estimate the contribution 
from a large-scale magnetic field. After subtracting out 
the estimated contribution from rotation, there is some 
residual signal in the even splitting coefficients. This may 
be explained by a magnetic field of approximately 20 kG 
strength located at r = O.96i?0 in 1996. Since the density 
here is of order 4 x 10~ 3 , and the downward velocity for 
the plumes is of order 5 x 10 4 cm s _1 (Nordlund 1999), 
the estimated dynamical pressure of the plumes, pv 2 , is 
equal to or larger then 10 7 dyne cm~ 2 . The size of pv 2 is 
comparable with the magnetic pressure, B 2 /8tt, for a field 
strength of 20 — 30 kG. Therefore, as a body, a magnetic 
flux tube experiences not only a buoyant force, but also 
a dynamic pressure of downflows. The total pushdown 
force of downward plumes equals to the dynamic pressure 
{pv 2 ) times the total cross section of the plumes (S). The 
dynamic acceleration g n can be expressed as follows: 

g n = 2sv 2 z ^/ita, (10) 

where s = 1 when v z > 0, s = —1 when v z < 0, 
£ = S/2aL < 1 is the fractional area of the downflows. The 
upward-moving plumes underneath the flux tube (s = 1) 
accelerate the magnetic buoyant diffusion of the tube, 
while the downward-moving plumes above the flux tube 
(s = — 1) go against the magnetic diffusion. These down- 
flows act as an anti-diffusion process. 

The magnetic field in a flux tube can be considered to 
be a "mean field" since the mean-field concept is relative 
in Mean Field Magnetohydrodynamics (e.g., Krause and 
Radlcr 1980). As a matter of fact, it depends on the spatio- 
temporal range over which the average is made. For the 
mean field in a flux tube, the spatio-temporal range is the 



volume and lifetime of a typical magnetic flux tube. With 
this in mind, the field in a flux tube obeys the classic dy- 
namo equation (Parker 1955), 



9B 

~dt 



= V x (U x B + £) - V x (77V x B). 



(11) 



Using this equation, we cannot study the formation of flux 
tubes from the large-scale magnetic field, but we can study 
if the flux tube fields can be amplified. 

We assume a pure rotation motion U = fi x r, where 
n is the angular velocity. This will cause the ^-effect. 
£ contains the a-effect, aB, and the magnetic diffusion 
processes, which are caused by the turbulent pressure, 
/3 t V x B, by the magnetic buoyant force, /3&V x B, and 
by the dynamic pressure of downward-moving plumes, 
/3„V x B.' Therefore, 

£ = aB-(p t +p b + /3„)V x B. (12) 

The magnetic diffusion caused by the microscopic gas col- 
lision motion, has been included in Eq. (11) in terms of the 
classic magnetic diffusivity 77V x B. For convection, a can 
be approximately expressed as (Krause & Radler 1980) 

a = \hT cor , (13) 



where h = v • (V x v) is the mean helicity over the corre- 
lation time T cor . 




Fig. 1. — Solar rotation rate as a function of radius and latitude 
(Antia ct al. 1998). 

In order to work out the various magnetic diffusitivites 
mentioned above, we observe that the diffusivity of a force 
is proportional to its impulse per unit mass. The impulse 
is defined as the force times its action time. The impulse 
per unit mass has the dimension of acceleration times time. 
The action time can be taken to be the propagation time 
of a magnetic perturbation across the diameter of the flux 
tube. Since the characteristic propagation velocity of mag- 
netic perturbations is the Alfven velocity Va, the action 
time is estimated to be i ac t = 2ci/Va- However, diffusivity 
has the dimension of acceleration times time times length. 
The charateristic length for a flux tube is its length L. We 
use this length L to meet the need of dimensionality for 
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the buoyant and the dynamic diffusivity. Consequently, 
we have 

fi b = g b (2a/V A )L = (4np)^ 2 2agLB / (8nP - B 2 ), (14) 

fin = g n (2a/V A )L = (Anp^Asv^L/irB. (15) 
The magnetic turbulent diffusivity fit = ^v' ^cor is known 
(Krause & Radler 1980) for an isotropic turbulence, where 
v' is the turbulent velocity. 
In summary, we define 

fi = 7] + fi t +fi b + fi n , (16) 
where fib and fid are given by (14) and (15), respectively. 
As a result, the hydromagnetic induction equation be- 
comes the dynamic equation of a magnetic flux tube, 
BTi 

— = V x (U x B + aB) - V x (fiV x B). (17) 
Since fi depends on B, this equation is nonlinear. 

3. LINEAR ANALYSIS 

The flux tube model equation obtained in the previous 
section is a nonlinear partial-differential vector equation. 
Since such equation is difficult to solve analytically, while 
a simple linear analysis may reveal some important char- 
acteristics about its solution, we analyze and solve the 
equation in the linear approximation in this paper. Unlike 
numerical simulations, the linear analysis cannot take into 
account full nonlinearity though it reveals some nonlinear 
effects such as critical magnetic fields of a flux tube. 

3.1. Linearization of equations 

Substituting Eq. (1) into Eq. (17), and linearizing the 
result with respect to the perturbation of B and A near B 
and A = 0, denoted by B' and A', we obtain two linear 
equations that govern B' and A', 

™-=aB' + fi V 2 A\ (18) 

^ = to™ - Cl™ - aV 2 A> + fi Q V 2 B> 
at ox oz 

+fi B B' - fi z ^- + fi x ^- (19) 
oz ox 

in the local Cartesian frame e z = e r , ej, = e^, e x = eg. 
We have defined, 
d 

Q x = —(rtosmO) = fl(cos9 + ilgsm9), 
Ox 

d 

tt z = —(rflsme) = Qsm9(l + fl r ), 



dz y 

fi B = fiW 2 B Q 



dfi dB dfi dB 



dz dz 
dB 



+ 



dx dx 



(20) 



Ox Ox 



fiz 



dfio 
dz 



■fi 



dB 

dz 



where fi = fi{B ,v z ), tog = d\nfl/d9, to r = d In to/d In r, 
and fi =_{dfi/dB)\ B=Bo = [fi, (B ) - fi 2 (B ,v z )]/B . The 
second line of Eq. (19) stands for the nonuniform effect. 
When this effect is ignored, Eqs. (18-19) are the same as 
the classic dynamo equations, except that fio depends on 
the magnetic field and the convection velocity field v z . In 
the classic case, fio is always positive, implying a genuine 
diffusion. In contrast, when the dynamic pressure of the 
downward moving plumes is larger than the buoyant pres- 
sure of the flux tube, fio can become negative. 



3.2. Dispersion relation 

In order to investigate the linear instability, we try the 
wavelike solutions of the form 



B' = B'o exp[i(cjt - k-x)], 
A' = A'o exp[z(o>i - k- x)], 



(21) 
(22) 



where u> is complex, while the other variables are real. 
Substituting Eqs. (21-22) into Eqs. (18-19), we obtain 



{fiok 2 + iuj)A' = aB , 



(23) 



(fiok 2 +iu-fi B + ifi k )B' = -i(k x Cl z - k z n x )A' ,(2A) 

where fik = fi x k x — fi z k z . The dispersion equation can be 
obtained by multiplying Eq. (23) and (24), 

(fiok 2 + iui)(fiok 2 + iuj - fi B + ifi k ) = -ia(k x Ct z - k z (l x ). 

(25) 

Solving this equation for to as a function of wavenumber 
vector k, we obtain the dispersion relation, 



u = u(k x ,k y ,k z ;v z ,B ), 



(26) 



where B and v z are parameters. 

Under the homogeneity approximation, we can neglect 
fi B and fik- In this case, Eq. (25) reduces to 

(fiok 2 + iuj) 2 = -ia(k x h z - k z Q, x ). (27) 

Rewriting —i as exp(i37r/2), when 

D = a(k x Cl z - k z Cl x )/2 > (28) 

is satisfied, we can solve this equation for u) easily as fol- 
lows: 

(29) 



w ± = ±\ r D + i(fik 2 Tv r D)- 



3.3. Period 

When D^O, the magnetic flux tube generation is pe- 
riodic. The period is 



T = 2n/VD. 



(30) 



It is well-known that a is positive (negative) for cyclonic 
convection in the northern (southern) hemisphere of the 
Sun. Therefore, the condition D > leads to 



(1 + tt r )k x > (cot 9 + Vlg)k z 



(31) 



When this condition is satisfied we have two magnetic flux 
tube waves propagating in the opposite dirrection. These 
two waves have the same period. 

In order to calculate T, we first determine k x and k z . 
The fundamental wavenumbers are, k x0 = V^ir/r in the 
x-axis, and k z0 — 2n/r in the z-axis at r since the max- 
imum wavelengths are A^o = V2r and A^o = r - The 
butterfly diagram implies that both branches share the 
whole length, so we assume k x — ±k x o for the equator- 
ward/poleward branch in the northern hemisphere of the 
Sun. The condition D > leads to kf = ±k z0 in the same 
hemisphere. 

Helioseismic data allow us to probe the rotation 
rate in the solar interior as a function of radius 
and latitude (Thompson et al. 1996; Antia et al. 1998; 
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Schou et al. 1998; Antia & Basu 2000). We use the rota- 
tion rate inferred from Global Oscillation Network Group 
(GONG) data by Antia, Basu, and Chitre (1998), as shown 
in Fig. 1. 

Table 1 shows T (unit: year) as a function of radius 
and latitude, calculated by using Eq. (30). This table 
shows a low latitude branch and a high latitude branch. 
Since the observed length of the solar cycles is equal to 
about 22 years, this table shows that the cycle should 
originate from the surface layer with r > O.95i?0, which 
is in agreement with the inferred depth from helioseismic 
data by Antia, Chitre and Thompson (2000), who found 
that the solar magnetic field peaks at r = 0.96Rq. The 
period also depends on latitude. If we use the observed 
period (22 years) as a criterion, this table indicates that 
the low- (high-) latitude branch originates from the lat- 
itude 35° (40°). Since k x > and k z > for the low- 
latitude branch, the magnetic flux tubes near the latitude 
35° and radius r — 0.96 Rq move upwards and equator- 
wards. Since k x < and k z < for the high-latitude 
branch, the magnetic flux tubes near 40° and 0.96i? Q move 
downwards and polewards. This scenario is in agreement 
with observations (Makarov & Sivaraman 1989). 

For different branches we have to use different a values 
to produce the desired period at the mid latitude. Since 
a is determined by the Coriolis force and the higher the 
Coriolis force the lower the latitude, it is reasonable to 
use a larger a for the low-latitude branch. Nevertheless, 
a jump of a from 90 cm s _1 at 35° to 5 cm s _1 at 40° is 
implied and needs to be explained. 

Table 1 also shows that solar magnetic activity is a 
multi-period phenomenon, and the longest period is 22 
years, occuring in the latitude range from 35° to 40°. 
The period decreases equatorwards, polewards, and down- 
wards. The period can be as small as 1-2 year near the 
base of the convection zone and near the equator. 

3.4. Critical magnetic field 

Eqs. (21-22) show that spontaneous generation of a mag- 
netic flux tube requires loi = ((3 k 2 T VD) < 0, which 
implies 

Po < ±v r D/k 2 . (32) 

For simplicity, we further neglect the classic collision and 
isotropic turbulent diffusivity. As a result, the unstable 
condition (32) becomes 



Pb + p n T v^D/k 2 < 0. 



(33) 



We assume Bo <C B g , where B g = (8ttP) ' is the mag- 
netic field strength when the magnetic pressure equals 
to the total pressure P. In this case inequality (33) is 
quadratic since Pb w {4irpy/ 2 2agLB/B 2 . Solving this 
quadratic inequality for B , we obtain 

B+ < B < B+ 

for the u + branch, where 

1/2 



(34) 



B ± = B g 



± 



DP 



8pa 2 L 2 g 2 k i 
DP 2s£v 2 



8pa 2 L 2 g 2 k 4 nag 



1/2- 



(35) 



is the critical magnetic field for the u + branch. When B± 
is complex, we set B± = 0. For this branch, we have both 
upper and lower critical magnetic field. The former is non- 
negative, while the latter can be negative. The fact that 
Bo represents the amplitude of the toroidal field requires 

(36) 



Bo >0. 
we reset BZ' 



Therefore, when B^ < we reset BZ = 0. When 
B~\_ = BZ = 0, this branch vanishes. When v z < 0, 
s = — 1. In this case BZ < 0. So we reset it to be zero. 
Similarly, solving inequality (33) for Bo, we obtain 



BZ < Bo < B 



+ • 



(37) 



for the lo branch, where 



b: 



b„ 



± 



DP 



1/2 



8pa 2 L 2 g 2 k 4 

2s^v 2 



DP 



8pa 2 L 2 g 2 k 4 nag 



1/2 



(38) 



When B± is complex, we set B± — 0. For this branch, 
we also have both upper and lower critical magnetic field. 
However, the latter is always non-positive, so we reset 
BZ = 0. The former can be negative, too. When BZ < 
we reset BZ = 0. When BZ = BZ = this branch van- 
ishes since B > is required. Only when v z < is BZ 
positive. 




Fig. 2. 

latitude. 



Critical magnetic field as a function of radius and 



From Eqs. (35) and (38) we can see that the critical mag- 
netic field of a magnetic flux tube depends on not only the 
velocity (including direction s and speed \v z \) of flows, the 
size (including radius a and length L) of the flux tube, 
and its location, but also the fractional area £. When the 
afi-effect is neglected (letting D — 0), the critical field 
is determined by competition between buoyant force and 
dynamic pressure of downflows, 

B% = B^-Zsivl/nag) 1 ' 2 . (39) 
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This is equivalent to the magnetic field by equating the 
buoyant and dynamic acceleration given in Eq. (7) and 
Eq. (10). Obviously, the critical field does not vanish only 
when there is a downflow so that s = — 1 and £ ^ 0. In 
fact, sunspots are often observed near the edges of gran- 
ulations where there exist downflows. The critical field is 
proportional to the speed of downflows, but inversely pro- 
portional to the square root of the radius of the flux tube. 
Fig. 2 shows as a function of solar radius and latitude 
with a = 2000 km, v z = -0.5 km s" 1 when r > 0.95i? Q , 
£ = 1, and L = 50a. The other parameters are the same 
as in Table 1. 

When the aO-effect is included, the critical magnetic 
field for the equatorward branch is not zero even if there 
is no downflow. In this case, Eq. (35) reduces to 



B„ 



DP 



2pa 2 L 2 g 2 k i 



1/2 



(40) 



Obviously, the thinner and shorter the flux tube the larger 
the critical magnetic field. The critical field without down- 
flows is weaker than the corresponding field with down- 
flows, but it does not vanish. The reason is that afi-effcct 
goes against the buoyant instability. In this case, the crit- 
ical field for another branch is zero. 



3.5. Growth rate 

From the dispersion relation given by Eq. (29), we find 
the growth rate for the u;* branch, 



-uf = -((3 b + f3 n )k z ±VD 
V2agLk 2 (B% - B ){B - B±) 
c B g \Bo\ : 



(41) 



where c = (P/p) 1 / 2 is the sound speed, and B± are given 
in Eqs. (35) and (38). Obviously, the growth rate is posi- 
tive when 



B ± <B < B±, 



(42) 



zero when B = B±, negative otherwise. This predicts 
that a magnetic flux tube experiences an ascending phase 
and a descending phase, implying that sunspots grow 
slowly, last for some time, and then disappear, since a 
sunspot appears when a flux tube crosses the surface of 
the Sun. This is in good agreement with sunspot observa- 
tions. 




Fig. 3. — Growth time of the critical magnetic field shown in 
Fig. 2 as a function of radius and latitude. 

Eq. (41) shows two wave branches, one that propagates 
equatorwards, and another that propagates polewards, can 
be excited. Fig. 3 shows the required time 

r B + 

r± = / [BY ± {B)]- 1 dB (43) 

for the tube field to grow from B\ = 10 _3 B^ to B+ (given 
by Fig. 2) as a function of solar radius and latitude. 
When v z = 0, Eq. (41) reduces to 



= ±v r D^/3 b k 2 . 



(44) 



This shows the magnetic field in a magnetic flux tube can 
grow since the afi-effect can go against the buoyant in- 
stability in the equatorward branch even if there is no 
downflow. Nevertheless, the growth rate is much smaller 
than that with downflows. 

3.6. Dynamo waves 

Solving Eq. (23) for A' , we can express the temporal 
evolution of the toroidal magnetic field B and the poloidal 
vector potential A in terms of the initial field Bo for each 
magnetic flux tube as follows: 

'27rf , 

=F k • x 



B* 



-Doe 



exp 



4ir 



exp 



T 



. 2nt , 7T 
z[— Tk-x-- 



(45) 
(46) 



The wavenumber vector k should satisfy Eq. (28). Substi- 
tuting Eqs. (45-46) into Eq. (1), we obtain, 



B± = B Q e r±t 



exp 



27rf , . 



{e 



cxT r 7r 

+ — (-k z e x +k x e z )exp -i-(l±2) 

47T L 4 



(47) 



where k x = y/2n/r, k z — 2tt jr and k y = are required by 
the condition D > 0, and ^ are the initial phase angles 
for these two branches. 



7 



This equation shows that the toroidal field leads the 
poloidal field by a phase of 37r/4 for the equatorward 
branch, but lags the poloidal field by a phase of 7r/4 for 
the poleward branch. The former is in agreement with the 
observations, which reveal that the poloidal component is 
almost in antiphase with the toroidal component. 

From Fig 3 it can be seen that the required time for 
the magnetic field of a typical magnetic flux tube to grow 
from a very low level to the critical field is very short 
in comparison with the cycle period when downflows are 
present. Using Eq. (43) it is easy to check that the decay 
time for such a magnetic flux tube is of the same order 
of magnitude as the growth time. Observations show that 
sunspots last, on average, a week or two. Therefore, our 
theory and observations are in agreement with each other. 
However, this timescale is much shorter than the cycle pe- 
riod. This suggests that magnetic flux tubes are only the 
basic elements of the magnetic flux tube wave described 
by Eq. (17). Therefore, we use the total magnetic flux 
$ = Nira 2 B rather than the magnetic field strength B 
for each magnetic flux tube as the wave variable of the 
magnetic flux tube wave, where N is the total number of 
magntic flux tubes present at the same time, a is the sta- 
tistically averaged radius of magnetic flux tubes, and B is 
the statistically averaged magnetic field strength of mag- 
netic flux tubes. Since a and B are mean quantities, only 
the number of flux tubes is variable, we finally find that 
the flux tube number varies with time and latitude 



N = N exp 



2irt 



T k • x + , 



(48) 



since x = r[l + tan(7r/4 — $)]/V2, where $ is the latitude. 
No is a constant in our single-tube theory. 

For the equatorward branch, we can investigate how its 
peak propagates with time by setting 



27rf/T - k x x - k z z + (f> + = (n - 1)tt 



(49) 



where n equals to the integer part of 1 + 2t/T. Since 
k x — \f2n/r and k z — 2ir/r, we obtain the migration 
equation of the wave peak 



4 



tan 



2t 
T 





407rV 


(5 


~ 180/. 



(50) 



where we have assumed $ + = 40° when t = 0. Obviously, 
the integer n plays a role of cycle number. 




40 

Time (years) 
Fig. 4. — Calculated Butterfly diagram. 

For the poleward branch, we can only see the reflected 
wave because it propagates into the interior of the sun. 
The reflection point is assumed to be located at that point 
where the gas density is one order of magnitude larger than 
that of the source region (r src = 0.96i?©). This leads to 
r r fl w 0.86i?©. The phase velocity of the wave in the radial 
direction approximately equals to 100 cm/s. Therefore, it 
takes about to = 5-5 years for the wave to travel from 
the source region (r src = 0.96-Rq) to the reflection point 
(r r fl = 0.86i? Q ), and then to the surface (r = Rq). As a 
result, the wave peak propagation equation for the pole- 
ward branch reads as follows: 

2iT(t + t ) /T + k x x + k z z + <jr = (n' - 1)tt, (51) 

where n' equals to the integer part of 1 + 2(t + to) /T. The 
solution of this equation is 

~2{t + t ) 



$ = — -Ran 



— (n' — 1) — tan 





45ttY 


(5 


~ 180/. 



(52) 

where we have assumed $ = 45° when t + to = 0. This 
implies a phase lag of ir/2 with respect to the equatorward 
branch, as shown in Fig. 4, which plots Eqs. (50) and (52). 

Since No is a constant in our single-tube theory, we pre- 
dict that the total number of flux tubes, on average, does 
not vary with time. This is in conflict with observations, 
which show that the sunspot number increases in the as- 
cending phase, and then decrease in the descending phase 
of the solar cycle. This problem may be solved by taking 
into account interactions between magnetic flux tubes. 

4. COMPARISON WITH OBSERVATIONS 

In this section we try to interpret the relevant observa- 
tional facts 1 — 10 listed in §1 in order, using the results 
obtained in the previous section. 

1. From Table 1 we can see that the cycle period 
equals to the observed 22 years when the dynamo 
operates above r = 0.95Rq, as inferred from 
helioseismology by Antia et al. (2000). Eqs. (30) 
and (28) show that the cycle period is determined 
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by the internal differential rotation and the 
cyclonic turbulent convection motions in the 
convection zone. 



2. Fig. 4 shows that the calculated butterfly diagram 
is in good agreement with the observed extended 
butterfly diagram by Makarov & Sivaraman 
(1989). 



What we observe at high latitudes are the reflected 
waves. They must be weaker than the original 
waves since the relfection is not necessarily 
complete. This naturally explains why the 
poleward branch is weaker than the equatorward 
branch. In other words, strong activity should be 
confined to low heliographic latitudes |$| < 35° of 
the equatorward branch, as observed. 



4. Phase dilemmas: 



(a) Eq. (47) shows that the toroidal field leads 
the poloidal field by a phase of 37r/4 for 
the equatorward branch. This is in good 
agreement with the observation: the poloidal 
field is almost in antiphase with the toroidal 
component; 

(b) The 7r/2 phase lag of the poleward branch 
with respect to the equatorward branch can 
be explained by the time lag of 5.5 years for 
the former to travel from the source region 
to the reflection point, and then from the 
reflection point to the surface. 



5. From Eqs. (35), (38), and Fig. 2, we can see that 
the allowed magnetic field (i.e., the critical field) in 
a concentrated flux tube (specified by cross section 
radius a and tube length L) is strong at the depth 
of r = O.96i?0. Fig 3 shows that the growth (or 
decay) time for a typical flux tube is about several 
weeks at the depth. 



6. We use the internal rotation rates obtained from 
helioseismology (Fig. 1). 
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Fig. 6. — Assumed solar internal magnetic field as a function of 
radius in 1996, as suggested by helioseismology (Antia et al. 2000). 

As elucidated by Li & Sofia (2001) and Sofia & Li (2001), 
the cyclic variations of global solar parameters (Fig. 5 ) can 
be reproduced by the magnetic field peaked at the depth 
of r = 0.96i? Q depicted in Fig. 6: 

8. The first and third panels of Fig. 5 show that the 
total solar luminosity varies in phase with the solar 
cycle. The amplitude equals about 0.1%. 



7. Fig. 2 shows that the magnetic fields peak at 
r = O.96i?0 when high-speed downflows are 
assumed to occur above r = O.95i?0. Though 
the critical field is stronger at the base of the 
convection zone, the growth time is longer than 
the cycle period. 



9. The first and fourth panels of Fig. 5 show that the 
solar effective temperature varies in phase with the 
solar cycle. The amplitude equals about 1.5°C. 

10. The first and fifth panels of Fig. 5 show that the 
solar radius varies in phase with the solar cycle. 
The amplitude equals about 20 mas. 
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If the magnetic field is placed at the base of the convection 
zone, we cannot find a strongth that simultaneously repro- 
duce the observed cyclic variations for these three global 
solar parameters. 

5. CONCLUSION 

We use 10 main observational facts pertinent to the solar 
cycle activity to filter the CZ, OL, IF and BL solar cycle 
models. We find that none of them are satisfactory. In par- 
ticular, recent helioseimic inversions (Antia et al. 2000) 
that determine the location of solar internal magnetic 
fields at r = 0.96i? Q , rule out the OL, IF, and BL mod- 
els, which place the solar magnetic fields at the base of 
the convection zone. This is one of the main motivations 
to reconsider the CZ dynamos. The second motivation 
comes from trying to reproduce the observed cyclic varia- 
tions of total solar irradiancc (Frohlich & Lean 1998), so- 
lar effective temperature (Gray & Livingston 1997), and 
solar radius (Emilio et al 2000). Extensive numerical ex- 
periments (Li & Sofia 2001) shows that only magnetic 
fields located in the convection zone can simultaneously 
produce the observed cyclic variations of these three 
global solar parameters. The third motivation originates 
from numerical simulations of the solar convection zone 
(Chan and Sofia 1989; Nordlund 1999) because they indi- 
cate a major presence of downward- moving plumes of high 
velocity. 

For the CZ dynamos, we have to face the magnetic buoy- 
ancy instability problem. The best way to do so is to in- 
corporate the buoyancy into the dynamo cuqations. This 



forces us to treat the mean field in a flux tube. To include 
the dynamic pressure of a flow, the flux tube should be 
considered to be a body. As a field, both the cyclonic con- 
vection (a-effect) and differential rotation (fi-effect) play 
a role. As a body, the tube experiences not only a buoyant 
force, but also a dynamic pressure of downflows above the 
tube. We show that these two dynamic effects can be in- 
corporated into the dynamo equations by adding to them 
two diffusion terms. 

We analyze and solve the extended dynamo equations 
in the linear approximation by using the observed solar in- 
ternal rotation rates and assuming a downflow suggested 
by numerical simulations of the solar convection zone. The 
results are in agreement with all the properties of the solar 
cycle listed in §1. 

We can not only interpret the observational facts of the 
22-year cycle, but also make interesting predictions. For 
example, Figs. 2 and 3 may imply longer activity periods 
of dacades to centuries in addition to the 22-year cycle, 
as suggested by the historical data. We'll explore them 
separately. 

We want to thank Dr. Sarbani Basu for her kindly pro- 
viding us with rotation rate data, Dr. Frank J. Robinson 
for his help, and Prof. Pierre Demarque for useful discus- 
sions. This work was supported in part by a grant from 
the National Aeronautics and Space Administration, and 
in part by Natural Science Foundation of China (project 
19675064). 
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Table 1 

The length of the solar cycle (unit: year) as a function of radius and latitude. 
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